Repulsive and attractive Casimir interactions in liquids 
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The Casimir interactions in the solid-liquid-solid systems as a function of separation distance 
have been studied by the Lifshitz theory. The dielectric permittivity functions for a wide range of 
materials are described by Drude, Drude-Lorentz and oscillator models. We find that the Casimir 
forces between gold and silica or MgO materials are both the repulsive and attractive. We also 
find the stable forms for the systems. Our studies would provide a good guidance for the future 
experimental studies on the dispersion interactions. 
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I. INTRODUCTION 

The dispersion interactions, the Casimir force, between 
neutral objects have brought attraction for many years. 
There are a lot of factors affecting on the value of force, 
such as geometry and material properties. Each of them 
gives rise to hot subjects of ongoing investigation. Some 
experiments have examined the influence of the dielec- 
tric properties of objects on the Casimir force [TH3J [6] . A 
number of settings used to study the interaction in terms 
of theory are ideal metals, real metals and semiconduc- 
tors [2 |H [6], metamaterials, and two objects placed in 
liquids [UEMl]. These studies have significantly advanced 
our understanding of the subtle effect of geometry and 
material on the Casimir-Lifshitz interactions, especially 
for designing nanodevices and nanotechnologies. 

In the Lifshitz theory, the dispersion interactions pri- 
marily depend on dielectric permittivity functions of ma- 
terials. Changing dielectric function alters the Casimir 
interactions. There are some ways to modify dielectric 
functions, including illuminating a light on the silicon 
[H] , which make drifting carriers on semiconductor ma- 
terials. In principle, there are some models to describe di- 
electric response functions of real materials, for example, 
plasma and Drude models for metals [51 [HI E] > Drudc- 
Lorentz and oscillator models for liquids [2J [TT] , oxides 
and others [IOlfT2] . Based on these models, the Casimir 
forces were obtained by numerical integrations and series 
expansion methods [T5] . 

It has been theoretically shown that the attractive 
Casimir interaction always occurs between two (non- 
magnetic) dielectric bodies related by reflection. There- 
fore, the repulsive force is a striking feature creating in- 
spiration for scientists to make accurate measurements 
of nano electromechanical machines where the repulsive 
force plays an important role and might resolve the stic- 
tion problems. The repulsive Casimir forces can be ob- 
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served in systems which have the presence of liquids [2] , 
metamaterials and metallic geometries [IB]. Recent ex- 
periments have pointed out that there a repulsive force 
exists between a gold sphere and a silica plate, separated 
by bromobenzen [2]. As a matter of fact, the repulsive 
Casimir forces between solids arise when the dielectric 
of material surfaces 1 and 2 and an intervening liquid 
obey the relation £i(i£) > £u qu id{iC) > £2^0 over a 
wide imaginary frequency range £. 

A previous theoretical 14 study has noticed that it 
is difficult to establish an equilibrium configuration of 
sytems in a vacuum medium. In the reference [15] . the 
authors showed that they were able to form some sta- 
ble configurations of Teflon-Si and Silica-Si immersed in 
ethanol. The equilibrium is explicitly explained by dis- 
persion properties. In the present work, our theoretical 
studies have shown that the equilibra can be obtained 
by placing Au-MgO, Silica-MgO and Au-Silica sytems in 
bromobenzen. 

In this paper, the Casimir-Lifshitz forces between ma- 
terial plate systems made in oxides and metals immersed 
in bromobenzen are calculated. The combination be- 
tween these results and the proximity force approxima- 
tion (PFA) method allow us to compute the Casimir in- 
teractions in different configurations. We find that the 
magnitude of the Casimir force between two dielectric 
bodies depends on the configuration and distance be- 
tween two bodies. The shape usually used in experiments 
is a combination of a sphere and a plate because one can 
avoid the problem of alignment and easily control the dis- 
tance between them. The energy interactions between a 
plate-plate system per unit area can be obtained by us- 
ing the relationship between the Casimir energy of two 
plannar objects and the dispersion force of a sphere-plate 
system. 

The rest of the paper is organized as follows: In Sec. 
II the theoretical formulations of Casimir-Lifshitz force 
interaction are introduced. In section III, the numerical 
results for the Casimir force between two bodies are pre- 
sented. Important conclusions and discussions are finally 
given in section IV. 



2 



II. LIFSHITZ THEORY FOR FORCE 
CALCULATIONS 

For the force calculations, we used Lifshitz theory with- 
out considering effect of temperature. The separations 
used here were less than 1 /xm, therefore thermal cor- 
rections at T = 300 A' are not significant. As previously 
noted in [2, 3, 20, 21 , the Lifshitz formula at zero temper- 
ature for the Casimir force acting between between two 
parallel flat bodies per unit area, separated by a distance 
d are given by 
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Here the reflection coefficients \ t TE and r^L te f° r 
two independent polarizations of the electromagnetic 
field (transverse magnetic and transverse electric fields) 
are 
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in which e^ p ' (if) and (if) are the dielectric permit- 
tivity and the magnetic permeability of the first body (p 
= 1) and the second body (p=3), respectively, e^ 2 ^ (cj) 
and ^ (if) are the dielectric function and the per- 
meability of a liquid filled between two bodies. Here, 
medium 2 selected is a bromobenzen so /i^ 2 ' (if) = 1. 
Moreover, in this paper, the non-magnetic materials used 
such as germanium, gold and oxides have also /j,^ (if) = 
1. k± magnitude of the wave vector component perpen- 
dicular on the plate, is frequency variable along the imag- 
inary axis (u = if). 

We recall that Lifshitz formula, routinely used to in- 
terpret current experiments, express the Casimir force 
between two parallel plates as an integral over imagi- 
nary frequencies if of a quantity involving the dielectric 
permittivities of the plates us = if. It is important to 
note that, in principle, recourse to imaginary frequen- 
cies is not mandatory because it is possible to rewrite 
Lifshitz formula in a mathematically equivalent form, in- 
volving an integral over the real frequency axis. In this 
case, however, the integrand becomes a rapidly oscillat- 
ing function of the frequency, which hampers any possi- 
bility of numerical evaluation. Another remarkable point 



is that occurrence of imaginary frequencies in the ex- 
pression of the Casimir force is a general feature of all 
recent formalisms hence extending Lifshitz theory to non- 
planar geometries jTTJ [T5] . The problem is that the elec- 
tric permittivity e(if) at imaginary frequencies cannot 
be measured directly by any experiment. The only way 
to determine it by means of dispersion relations, which 
allow the expression of e(if) in terms of the observable 
real-frequency electric permittivity e(if). In the stan- 
dard works on the Casimir effect, e(if) is expressed with 
the Kramers-Kronig relation in terms of an integral of 
a quantity involving the imaginary part of the electric 
permittivity |13) 
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where Ime(w) is calculated using the tabulated optical 
data for the complex index of refraction. 

The well-known dielectric function described for gold 
is the Drude model [13] 
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where u p — 9.0 eV, 7 = 0.035 eV are the plasma fre- 
quency and the relaxation parameter of Au, respectively. 

The imaginary part of the resulting dielectric function 
at 6 and 295 K of pure MgO are shown in [12 . The 
optical features have been fitted to a classical oscillator 
model using the complex dielectric function 
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where £00 is a high-frequency contribution, and UTo.ji ^1% 
and io V j are the frequency, full width and effective plasma 
frequency of the jth vibration. The values of these pa- 
rameters can be found in |12j . Of course with such simple 
model for the permittivity of MgO, there is no need to 
use dispersion relations to obtain the expression of e(if ), 
for this can be simply done by the substitution uj — > if 
in the r.h.s of Eq.([8|) [21 
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In the case of bromobenzen and silica, it has recently 
been used for measurement of repulsive forces between 
gold and silica surfaces. Extremely weak repulsion was 
measured, indicating that the dielectric functions of bro- 
mobenzen and silica are very similar in magnitude. In 
fact, oscillator models are constructed to represent the 
dielectric function at imaginary frequencies. The form of 
the oscillator model is given by 
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where the coefficients C, are the oscillator's strengths cor- 
responding to (resonance) frequencies w, [U El 123 • The 
dielectric data was fitted in a wide frequency range [2J. 
They are more accurate in comparison with other simple 
oscillator models. Moreover, many older references used 
limited dielectric data, so the oscillator models with sec- 
ond or third order may lead to the difference in Casimir 
force calculations. The parameters we used in the present 
paper for bromobenzen and also silica come from [5] . 

III. NUMERICAL RESULTS AND 
DISCUSSIONS 

The Casimir attractive force usually occurs in experi- 
ments and theoretical calculations. When bromobenzen 
is filled in the gap between two bodies, the Casimir force 
is attractive if the dielectric functions do not satisfy one 
condition £i(i£) > £u q uid(i£,) > £2{i£,) for all frequencies 
£. Therefore, by describing Fig. [T] as the dielectric re- 
sponse function as a function of the frequency gives us 
some predictions of repulsive and attractive forces. 
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FIG. 1: (Color online) The dielectric function of various ma- 
terials plotted at imaginary frequencies £. 

This graph shows that £Au(iO > SMgo(i£) > 
£uguid(iO and e Mg o{iCl > £Au(iQ > EiiquidiK) at 
£ < 6.5 eV, thus the interactions between Au and MgO 
body immersed in bromobenzen liquid and in vacuum 
are attractive in this range. In the range of £ > 6.5 
eV, e Mg o(iO > £liquid(i£) > £A«(i£)> it; causes the re- 
pulsive interaction. Similarly, in the gold-bromobenzen- 
silica system, at extremely small frequencies , the forces 
are attractive. In the larger frequency region, the Casimir 
forces are repulsive. Besides, the similar explainations 
are applied to understand the interaction in the MgO- 
bromobenzen-Au system. The numerical calculations of 
the normalized Casimir force are provided in Fig. [2] 

In the MgO-bromobenzen-Silica system, it can be 
clearly seen that there are two positions in each curve 
where the Casimir force is equal to zero. The first 
points are corresponding to unstable equilibria dus ~ 13 
nm because the interaction force changes from the at- 




FIG. 2: (Color online) Relative Casimir force between two 
semi-infinite plates normalized by the perfect metal force 
F (d) = — Tv 2 hc/240d 4 . The liquid used in this calculation 
is bromobenzen. 

tactive force to the repulsive force, the second point 
da ~ 110 nm is a stable position. There is only one 
position in the Au-bromobenzen-Silica system and the 
Au-bromobenzen-MgO sytem, The interaction forces dis- 
appear at d s ~ 275 nm and d s rs 5.5 nm, stable posi- 
tion of each sytem, respectively. 




FIG. 3: (Color online) Schematic picture of the setting consid- 
ered in our calculations. A sphere is located in bromobenzen 
at a distance d away from a material plate. 

In order to consider the Casimir interactions between 
a spherical body and a plate at a distance of close ap- 
proach d at a temperature T = 300 K, it is very useful to 
utilize the PFA method to calculate. Experimental re- 
sults for the Casimir force in the plane-sphere geometry 
are usually compared with PFA-based theoretical mod- 
els. The spherical surface is assumed to be nearly flat 
over the scale of d. Although the Casimir force is not 
additive, PFA is often expected to provide an accurate 
description when R 3> d. Here, the radius of Au sphere 
that is used in configurations is R = 40 ^m in order to 
calculate Casimir interaction by the proximity force ap- 
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proximation (PFA) method because the ratio of d to R 
is small enough to PFA results becoming enormously ac- 
curate. It can be described by Fig. [3] In this approach, 
the surfaces of the bodies are treated as a superposition 
of infinitesimal parallel plates [22] , 
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here F pp is the Casimir force for two parallel plates of 
unit area. 

When using the PFA method, one important point is 
that the interactions between a gold sphere or a mag- 
nesium oxide sphere and a silica plate are equal to the 
interactions between a magnesium oxide plate, which has 
the same radius, and a gold plate or a silica plate. There 
is no difference in calculations and results as well because 
the PFA method does not consider a structure of bodies 
when their shape is modified or is spherical or cylindrical 
shape. The equivalent situations occur in other materi- 
als. The resulting Casimir forces are shown in Fig. [4j 
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FIG. 4: (Color online) The Casimir forces of various sphere- 
bromobenzen-plate systems are estimated as a function of sep- 
aration, described in the text with the spherical redius R = 40 
fim. 

In the reference [2], the authors experimentally mea- 
sured and theoretically calculated the Casimir interac- 
tion between a gold sphere and a silica plate immersed 
in bromobenzen in the range from 20 nm to 60 nm. Our 
results in this range are the same for this range. But 
when we extend the considered range of distance, the 
attractive-repulsive transition occurs at approximately 
190 nm. This position makes this system stable. Another 
consequence of Fig. [4] demonstrates that stable position 
of the Au-bromobenzen-MgO sytem moves to 3.5 nm to 
balance between the attractive and repulsive forces. It 
can be explained that increasing the separation distance 



of infinitesimal parallel plates causes the fast reduction 
of the dispersion interaction. At the same minimal sepa- 
ration distance d, the attractive force acting on a sphere 
is less than that of a plate in the same effective area. Fi- 
nally, in the system of a MgO sphere and a silica plate 
embeded in bromobenzen, there is only one presence of 
non-interaction posion at nearly 10 nm. It is an unstable 
position. 
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FIG. 5: (Color online) The Casimir energy is calculated as a 
function of separation for different materials . 



In addition, PFA formula and Eq.(ll) can allow us to 



estimate the Casimir energy per unit area between two 
plate bodies illustrated in Fig. [5] The Casimir energy is 
approximated by [22] 



F FFA (d) = 2irRE(d), 



(12) 



where E{d) is the Casimir energy per a unit area for 
planar bodies. 

We have also applied the PFA method to calculate the 
Casimir force in sphere-sphere systems, we continue to 
calculate by the PFA method. The formula for this cal- 
culation is given 

F FFA (d) = 2tt J rdr x 
o 



F pp (d + Rx - ^R\ - r 2 + R 2 - \jR\ - r 2 ), (13) 

where the radii of two spherical objects are R\ and R 2 , 
respectively. It is assumed that R 2 < R\. In this study, 
we consider R\ — 40 /im and the case of R 2 = Ri, Ri = 
2R 2 and R 1 = 2R 2 . 

Here, having calculated F FFA (d) in a sphere-sphere 
system using Eq.(fl3l and F F p FA (d) in a sphere-plate 
system using Eq.(fTT|). These results obtained show 
that, when increasing d, the ratio F FFA (d) / F F p FA {d) 
does not depend on the distance d. It is a constant 
with its magnitude as a function of the radius of two 
spheres, F FFA (d)/F FFA {d) = 1/2 when R ± = R 2 and 
Fss FA {d)/F F p FA {d) = 1/3 when R x = 2R 2 . Generalizing 
this ratio, if i?i = nR 2 , we have F FFA (d)/F FFA (d) = 
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l/(n + 1). This character is likely to be explained by 
the results in [22]. When the second sphere is extremely 
small in comparison with the first one, the interaction 
force goes to zero. In this case, the Lifshitz formula used 
to calculate the Casimir force should be transformed to 
the Casimir-Polder formula describing the interaction be- 
tween an atom and a microscopic object. Moreover, the 
PFA method is not accurate because this approach is use- 
ful if the size of objects is much larger than the separation 
between them. On the other hand, we have 
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where i?i = nR 2 . If Ff p FA = 0, F FFA must be zero. 
Therefore, the unstable and stable positions are constant 
and unchanged when a radius of a second sphere varies. 

One demonstrated that the Casimir force between two 
objects embedded in liquids can be derived from the well- 
known Lifshitz formula at least if the object is not made 
of nonabsorbing materials [5] . That explains why the Lif- 
shitz expressions is used in order to calculate the Casimir 
force and compare with experimental data. Nothing 
changes in the dielectric functions of bodies immersed 
in liquids. On the other hand, several experiments ver- 
ified that when metals are placed in liquids, there is a 
variation of Drude parameters in the metal |23j . The 
discrepacy of the interaction between "dry" and "wet" 
can reach to 15 % in this case. But they measured the 
Casimir force between two metal plates and got the error. 
Besides, maybe the dielectric of liquids and low index ma- 
terials play much more important role in Casimir force. 
In the reference [5] , the change of Drude parameters are 
not taken into account but the theoretical calculations 
are close to the experimental data curves when we have 
liquids and the low index materials. 



IV. CONCLUSIONS 



In this work, we have extended the Lifshitz theory 
to calculate the Casimir force. Liquid, silica and mag- 
nesium oxide are represented by oscillator models. Al- 
though further studies are required to determine the re- 
pulsive Casimir force accurately, our results show that 
MgO and silica is a good candidate for the demonstration 
of quantum lcvitation. The contribution of bromobenzen 
is important because it is an important factor making the 
purely repulsive force or the repulsive-attractive transi- 
tion. After calculating the Casimir force between two 
bodies per unit area and associating proximity force ap- 
proximation method, it is easy to compute the interaction 
of different material plates with a material sphere. Based 
on the formula of the Casimir force between a sphere and 
a plate, it is convenient to estimate the free energy in- 
teraction of bodies. The result is a prediction for further 
experimental studies. 
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